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Abstract. Progressive image transmission (PIT) is often used to reduce the transmission time 
of an image telebrowsing system. A side effect of the PIT is the increase of computational 
complexity at the viewer’s site. This effect is more serious in transform domain techniques 
than P in other techniques. Recent attempts to reduce the side effect are ; futile as they create 
another side effect, namely, the discontinuous and unpleasant image build-up. Based on a 
practical assumption that image blocks to be inverse transformed are generally sparse, t 
paper presents a method to minimize both side effects simultaneously. 


1. Introduction 

One important evaluation criterion for a telebrowsing system is the response time which i^ 
the time elapsed from the moment a retrieval request is issued until the desired information 
is actually displayed on the monitor [1], The response time can roughly be divided into three 
major parts. The first part is the searching time for the system to locate the desired 
information. The second part is the transmission time to send the information trough a 
channel. The third part is the display time for the information to be displayed on the 
monitor. The early studies of the telebrowsing systems were concentrated on the efficient 
retrieval of pure text information [2,3]* In this case, the searching time is the on y major 
concern However, for modern telebrowsing systems where multimedia information, 
including text, audio, image, and video, is considered, the transmission time and the display 
time become a significant part of the response time because of huge amount of data 
involved in still images and video (a sequence of images). 

To reduce the transmission time of an image telebrowsing system, a well known scheme 
called progressive image transmission (PIT) is often used. PIT allows an approxima e 
reconstruction of an image whose fidelity is built up gradually until the viewer decides either 
to abort the transmission sequence or to allow further reconstruction. This scheme increases 
the effective compression ratio because usually only a small part of the compressed data 
needs to be sent for browsing purpose. 

With PIT techniques, the transmission time can be greatly reduced. However, it also creates 
a side effect, that is, it increases the processing time at the viewer’s site because an inverse 
PIT process is required. Since the major task of the inverse PIT process is the image 
decompression given part of the compressed data, the research is aimed at the development 
of fast image decompression schemes for the inverse PIT process. 

The rest of the paper is organized as follows. First, the PIT schemes and their computational 
complexities are briefly addressed. Then, the drawbacks of recent attempts to reduce the 
computational complexities are discussed. Next, the demonstration of a new approach is 
given. Finally, a performance comparison between the new approach and the recent ones 

is made. 
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2. PIT Schemes and Their Computational Complexities 

There are many PIT schemes. According to Tzou’s classification, they are divided into three 
major categories, namely, spatial domain, transform domain, and pyramid-structured, based 
on where the progression takes place [4], Each category can be further divided into several 
classes of techniques. The classification is shown in Figure 1. Note that not all of the PIT 
schemes will produce a considerable amount of computational overhead in the inverse PIT 
process. For instance, the spatial domain schemes only require a very low computational 
effort in the inverse PIT process. In pyramid-structured PIT schemes, only successively 
filtered pyramid techniques require high computational complexity in the inverse PIT 
process. Even for the successively filtered pyramid techniques, however, the complexity to 
process the first few levels of a pyramid from the top remains low. From a practical point 
of view, the processing of the first few levels of the pyramid may suffice the purpose of 
image browsing. On the other hand, transform domain techniques usually take considerable 
amount of computation in the inverse PIT process, since the inverse transforms have to be 
carried out with about the same computational effort for every stage of image 
reconstruction. 6 


Spatial Domain Techniques 
Bit-plane method 
Tree-searched vector quantization 
__ Progressively quantized DPCM 
Transform Domain Techniques 

Scanning pattern techniques 
Transform domain multistage quantization 
Bit-slicing method 

Pyramid-structured progressive transmission 
Tree-structured pyramid 
Binary tree 
Quadtree 

Successively filtered pyramid 

Figure 1. Tzou’s classification of PIT schemes 

In transform domain PIT schemes, the transform coefficients are first quantized and then 
divided into segments. Only one segment of quantized coefficients is sent for one stage of 
image reconstruction. The only differences among all transform domain techniques are the 
ways to determine the segments and the order in which they are sent. One common feature 
among them is that the transform coefficients are only "partially" encoded where in a non- 
PIT or sequential scheme they are said to be "fully" encoded. 

For a transform domain non-PIT scheme, one MxN inverse transform for an M x N image 
block is needed. However, for a transform domain PIT scheme, r times of M x N inverse 
transform are needed for the image block, where r is the number of stages of image 
reconstruction. The lower bound for r is 1 but its upper bound depends on the image, the 
viewer, and the PIT scheme. Therefore, the computation load for inverse transform’ is r 
times heavier in PIT schemes than in non-PIT schemes. 

One transform domain scheme using discrete cosine transform (DCT) receives great 
attention, since the DCT has the energy packing capabilities and also approaches the 
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statistically optimal transform (i.e. Karhunen-Loeve transform) in decollating a signal 
ewerned by°a Markov process [5]. In addition, it is part of the recently approved JPEG 
standard [6 Jl. The JPEG standard has brought a tremendous impact on the image-coded 
related industry However, as far as implementation of the standard is concerned, the 
standard provides only a guideline. How to implement the standard 
application still relies on the ingenuity of designers For example, JPEG has _chose 
specify neither a unique forward DCT (FDCT) algorithm or a unique inverse DCT (IDCT) 
in its recommendation. This is because research in fast DCT algorithms is ongoing a 
single algorithm is optimal for all implementations [7]. Forthe application of inverse PIT, 
we vldll show that traditional fast two dimensional (2-D) IDCT algorithms can be accelerated 
to reduce the processing time at the viewer’s site. 


3. Previous Approaches and Their Drawbacks 

To relieve the computation burden of IDCTs in inverse PIT, the following approaches have 
been used. 

Approach 1: Use traditional fast algorithms for IDCT. The computational complexity is 
reduced from 0(N 4 ) by the definition of IDCT to 0(NTog,N) by traditional fast algorithm^ 
where N x N is the block size. There are many fast algorithms available for roCT. For an 
8x8 IDCT, one of the best algorithms reported so far takes 96 multiplications and 6 

additions [8]. 

Approach 2: Use a fast progressive reconstruction method. It is a combination of a special 
scheme and the use of approach 1. This approach was first proposed by Takikawa to 
perform fast progressive reconstruction for discrete Fourier transformed and Wa s 
Hadamard transformed images [9], Later, Miran and Rao followed ^similar denvahon 
by Takikawa and developed a fast progressive reconstruction tor DCT image: s [10]. Ihe 
basic idea of approach 2 is to decompose the N x N transformed block into lofcN + 1 
sparse matrices, each of which can be inverse transformed by 1 x 1, 2 x 2, 4 x 4,..., ana in 
N fast inverse transform algorithms. 

Approach 2 has some advantages over approach 1. First, the computational complexity is 
lower. For example, consider a 4-stage image reconstruction and an 8 x 8 image block_ 
Approach 1 takes four 8x8 IDCTs while approach 2 requires only one 1 x 1 IDCT one 2 
x 2 IDCT one 4 x 4 IDCT, and one 8 x 8 IDCT. The computational saving is obvious. 
Secondly, the delay time is reduced. The delay time is the time to wait for all the elements 
in a transformed block before an inverse fast transform can be performed. 

However, approach 2 has a serious problem, that is, it has a poor and discontinuous image 
build up. The reason is that the order in which the sparse matrices are formed and sent is 
not in the order of visual significance. In general, a DCT coefficient with higher variance 
(or energy) tends to be more visually significant than that with lower variance. It is we 1 
known that the DCT coefficient variances are highly correlated along the zig-zag scan [11 j. 
Approach 2 has a fixed transmission pattern that does not even close to the zig-zag scan 
This problem has been confirmed experimentally by Miran and Rao [10]. They ascribed the 
drawback to not having low frequency terms immediately adjacent to DC components in the 
intermediate stages of reconstruction. Another drawback of approach 2 is that it s i 
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requires all elements of the sparse matrices to start computing the inverse transform. Thus, 
the delay time is reduced but not eliminated. 

One more drawback for both approaches 1 and 2 is the computational redundancy of 
traditional fast algorithms in inverse PIT. If IDCT is used in image decompression, its input 
block contains only a few nonzero coefficients. In addition, if a PIT scheme is used, the 
input matrix to IDCT contains even fewer nonzero elements. To visualize the redundancy, 
consider the signal flowgraph of a fast IDCT algorithm. Since a zero presented at an input 
node contributes nothing to the output, the paths between a zero input node and output 
nodes are trivial or redundant. 

To get a picture on how many spatial frequencies retained on the average after the 
quantization, many 512 x 512 8-bit greyscale images and RGB components of color images 
were tested. In the test, JPEG’s coding scheme, including a recommended quantization 
table, was used. Part of the' test result is presented in column 3 of Table 1. It is shown, even 
for a very busy image such as baboon image, no more than a quarter of quantized 
coefficients are nonzero. Even with so many spatial frequencies set to zero, the 
decompressed images and their originals are perceptually indistinguishable. Next, consider 
the case when a small visible image degradation is allowed. To produce a small image 
degradation, the same test was repeated except that the round-off operation in JPEG’s 
scheme was replaced by truncation. The new numbers are shown in column 4. The 
decompressed images have only minor degradation, for it does not diminish our capability 
to recognize meaningful objects in the images. In fact, the image quality is good enough to 
be the last stage of PIT. The test shows that the number of nonzero quantized DCT 
coefficients decreases sharply at the minor expense of image quality. In the inverse PIT 
process, the average number of nonzero elements in an 8 x 8 matrix does not need to be 
higher than that in column 4. 


Table 1. Average # of Nonzero Quantized DCT Coefficients in an 8 x 8 Block 


Images 

Image Activity 

Round-off 

Truncation 

Lena 

Low 

6.13 

4.01 

Boat 

Medium 

9.20 

6.00 

Baboon 

High 

15.50 

9.80 


How much of the matrix must be zero for it to be considered sparse depends on the 
applications. Generally, a matrix is called sparse if there is an advantage in exploiting its 
zeros [12], It is well known that exploiting the sparsity can lead to enormous computational 
savings in many applications such as solving simultaneous equations with Gaussian 
elimination method. Inspired by this fact, it is curious to see if the sparsity of the input 
matrix can also be exploited to compute IDCT efficiently in the environment of inverse PIT. 
Since the characteristic of an input image block to IDCT is generally not considered in 
traditional algorithms, a nonconventional approach must be adopted to exploit the sparsity 
of the input matrix. The proposed approach will be presented in the following manner. First, 
we describe the goal to be accomplished by the approach. Then, the rationale of the 
approach is discussed. Next, based on the rationale, two methods are presented - one is too 
slow to be useful, the other is its fast version. The fast version is shown to be good enough 
for the practical use. 
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4. The Proposed Approach 


In the inverse PIT process, computation burden of IDCT and computation redundancy 
associated with traditional algorithms are two major problems. The inherent drawbacks in 
Takikawa’s or Miran and Rao’s approach present another problem in the inverse PIT 
process. In view of all these problems, our approach should meet the following goals. First, 
it must be fast and efficient. Second, it must allow a scanning pattern that can conform to 
the visual significance. Finally, it must have practically no delay time. 

For the ease of discussion, several terms are defined first. A target matrix is an image block 
consisting of the quantized DCT coefficients that are partially encoded for PIT. Performing 
an IDCT on a target matrix results in a matrix called goal matrix. The result of processing 
one nonzero element in the target matrix is called the partial contribution to the goal 
matrix. Throughout this paper, the partial contribution is treated as a matrix or all its 
elements depending on the context. 

Based on the definition of 2-D IDCT, only nonzero elements in the target matrix can 
contribute to the goal matrix. In fact, the value of each nonzero element can affect the 
values of all elements in the goal matrix. The idea of our approach is to completely ignore 
the zero elements in the target matrix and process each nonzero element separately and 
efficiently. The goal matrix is then updated periodically by adding the partial contribution. 
Therefore, the computation of IDCT is divided into two tasks, i.e., the computation of 
partial contribution and the update of the goal matrix. The idea adapts particularly well to 
the scheme where DCT coefficients are run-length coded (such as JPEG s). 


The definition of 2-D IDCT is 

„ M-l JV-1 


Sx> = 


E E c(«)c(v)F uv cos 


\j MN K -0 v=o 


(2x+ 1 )kti Y f (2y+ l)v7t 
2 M j { 2N 


where x = 0, ..., M-l, y = 0, ..., N-l, and 

c(k)=— if k= 0, 

v/2 

= 1 otherwise. 


( 1 ) 


The coefficient in front of the double summation of equation 1 is only a scale factor which 
requires essentially no computation (except a register shift operation) in practical 
applications where M = N and M = 4, 8, or 16 are often used. Thus, it is usually neglected 
when comparing the computational complexity among the fast algorithms of IDCT. By 
taking the scale factor out, equation 1 becomes 


M-l N-l 

4= E E c(k)c(v)F uv cos| 

u=0 v=0 


< 2 ** 1 ) “* W (2j '* 1)v ’ , 1 

2M J { 2N ) 


( 2 ) 


where 


f 

Jxy 2 &x y 


If [f ] W is defined as the partial contribution to the goal matrix [f x7 ] due to F uv alone, then 
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(3) 


ty, t = c(uMv)F„cos| . )cos( j 

The partial contribution can be obtained by the use of definition in equation 3 directly. 
Assume the values of cosine functions for different combinations of x and u are 
precalculated and stored as a table. The table can also be used as the values of cosine 
functions for different combinations of y and v with x and u replaced by y and v, 
respectively. Let Q be the number of multiplications required to find the partial contribution 
due to F uv . Then, Q=2MN if both u and v are not zero, Q=3MN if u = 0 or v = 0 but not 
both, and 0 = 0 if u and v are both zero. For an M x N target matrix with n (> 1) nonzero 
F uv , where n < < MN, the number of multiplications to get the goal matrix is from 2(n- 
1)MN to 3nMN. With this naive approach, no addition but 128 to 384 multiplications are 
required if M = N = 8 and n = 2. This is not good enough, since an 8 x 8 fast IDCT can take 
as low as 96 multiplications [8]. Therefore, a better way to compute the partial contribution 
is needed. 


Equations 2 and 3 are equivalent if only one term in the double summation of equation 2 
is nonzero. So the traditional fast algorithms for equation 2 can be applied to equation 3 
as well. However, the direct use of them to compute the partial contribution is not desirable 
since they contain high computational redundancy. We found that with a systematic 
reduction rule for the signal flowgraphs of traditional fast algorithms, a much faster way to 
compute the partial contribution than the naive approach is possible. The rule is based on 
the two attributes associated with the partial contribution, which we call the mirror effect 
and the reducible property. 


From equation 3, it can be readily shown that 


foyluv ~ ("1) [^M-l-x’.yluv 
[fxyOuv - (-l) V [f,N-l-/]uv 

Muv = (-l) u+v [fM-i-x’,N.,.y] 


•y Juv 


(4) 


where x’ and y’ are particular values of x and y, respectively. Equation 4 indicates that the 
partial contribution exhibits high degree of symmetry or mirror effect. Note that only 
possible sign changes are involved in equation 4 and practically no addition or multiplication 
is required. The significance of this result is that only a quarter of the partial contribution 
needs to be determined through additions and multiplications. The rest of them can be 
determined by simple copy operations and possible sign changes. 


The reducible property can be stated as follows. The M x N partial contribution due to F uv 
is equivalent to that of (M/m) x (N/n) partial contribution due to F u/mv , n , where cd(M,u) 
= m, cd(N,v) = n, and cd(a,b) is a common divisor between non-negative integers a and b. 
This statement can be proved easily by noting that 


M, N [g u ,v M/m,N/n[^xy]u/m ,v/n (5) 

where w N [f xy ] uv is the [y uv defined in equation 3. If cd(M,u) = l only, it is said to be an 
irreducible partial contribution in row. Similarly, if cd(N,v) = l only, it is said to be 
irreducible in column. If for some m>l or n>l, the partial contribution is said to be 
reducible. Note that the reducible property is separable, i.e, the reduction in row size and 
column size can be processed separately. The largest reachable reduction for F uv happens 
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when gcd(M,u) = m and gcd(N,v) = n, where m>l and n>l, and gcd(a,b) is the greatest 
common divisor between non-negative integers a and b. For example, if M = N = 8 ,u — 6 ,and 
v= 4, the 8 x 8 partial contribution due to F M is equivalent to 4 x 2 partial contribution due 
to F 31 , since cd(8,6) can be 2 and cd(8,4) can be 4. It is also a maximum reducible case for 
F m since gcd(8,6) = 2 and gcd(8,4) = 4. Note that gcd(a,0) = a. Therefore, if u and v are both 
zero, M x N partial contribution due to F uv is simply a 1 x 1 partial contribution due to F 00 , 
which is always Fqo/ 2 no matter what the values of M and N are. 

Combining the mirror effect and the reducible property can lead to a great saving in 
computation of partial contribution. Consider the following example: We want to compute 
the 8x8 partial contribution due to F^. Since gcd(8,4)=4, it can be reduced to the 2x2 
partial contribution due to F n . By using the mirror effect, only 1 x 1 of the 2 x 2 partial 
contribution needs to be determined explicitly, which is f 00 =F 44 /2. It can then be expanded 
to 2 x 2 partial contribution by the use of mirror effect: 

*00 *01 

*io *n 

where f 01 = -f 00 , f ip = -f 00 , and f n = f oo- Similarly, by using the mirror effect for F u , we expand 
the partial contribution to 4 x 4 : 


*00 *01 *02 *03 

*10 *11 *12 *13 

*20 *21 *22 *23 

*30 *31 *32 *33 

where f 07 = -f 01 , f ()3 = "*00’ *12~'*11> *13 = ‘*10’ *20 = ‘*10’ *21 = ~*11» *30~"*00’ * 31 - "* 01 > *22“ *11* *23“ *10’ 
f 32 = f 01 , f 33 = f^. Since cd(8,4) can be 2, the 8 x 8 partial contribution due to F 44 is equivalent 

to" the 4 x 4 partial contribution due to F 22 . By using the mirror effect for F 22 , we can expand 
the partial contribution of 4 x 4 to the desired result of 8 x 8. Note that no multiplication 
is required to determine the 64 elements of partial contribution due to F^. 

The basic principle to reduce the signal flowgraph of a traditional algorithm is by retaining 
only the nontrivial paths. This concept is demonstrated by an example. Consider the row- 
column or indirect approach of a fast 2-D IDCT for a 4 x 4 target matrix. Chen’s algorithm 
is chosen here because it is simple and well recognized [13]. Normally, 8 4-point 1-D IDCTs 
are needed to accomplish the task (with very complicated data reordering, 4 4-point IDCTs 
are enough [8]). However, in our case at most 3 4-point IDCTs are necessary (1 along the 
rows (or columns) of the target matrix to get an intermediate matrix and 2 along the 
columns (or rows) of the intermediate matrix to get a 2 x 2 submatrix of the partial 
contribution). The other three 2x2 submatrices can be derived automatically by the use of 
mirror effect. Furthermore, each 4-point IDCT can be done efficiently since only one input 
data out of 4 is nonzero. Consider the signal flowgraph for a 4-point IDCT shown in Figure 
2(a). The outputs of the 4-point IDCT (denoted by f 0 , f„ f 2 , and f 3 ) can be treated as linear 
combinations of the 4 inputs (denoted by F 0 , F„ F 2 , and F 3 ). Since only one of the inputs 
is nonzero. Figures 2(a) and 2(b) are functionally equivalent. The signal flowgraph in Figure 
2(b) can be further simplified by retaining only two of the four outputs (f 0 and fj) as shown 
in Figure 2(c) because the other two outputs can be derived by the use of mirror effect. 
Since the reducible property is separable, it can be used here to further reduce some of the 
subgraphs in Figure 2(c). Specifically, the subgraphs with input F 0 and F 2 are reducible. The 
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final irreducible subgraphs are shown in Figure 2(d). For convenience, the subgraphs shown 
in Figure 2(d) are said to be in their primitive forms. In other words, they can not be 
reduced or simplified any more. The above procedure can be extended easily to 8-point or 
higher order cases. 

The primitive subgraph with input F u and the one with input F v can be cascaded as a signal 
flowgraph to compute part of the partial contribution due to F uv . The connection rule is: at 
each output of the first subgraph, the second subgraph is cascaded. Which subgraph should 
be the first is immaterial as far as the result is concerned. However, the computational 
complexity may be different. 

The complexity of computing part of the partial contribution due to F uv can be examined 
by checking the primitive subgraphs with input F u and F y . The two primitive subgraphs are 
cascaded in the way described earlier. If the first subgraph takes P multiplications and the 
other requires Q multiplications, then the total number of multiplications required to obtain 
part of the partial contribution would be P + PQ multiplications. Alternatively, P and Q are 
also the number of output nodes for one subgraph and another, respectively. So P and Q 
can be obtained by counting the number of output nodes of the irreducible subgraphs. Since 
P+PQ = P(l + Q), a fast way to tell the required number of multiplications is to take the 
product of P and Q+l. P + PQ multiplications will also be the complexity to compute the 
full size partial contribution since no addition operations are involved and the expansion of 
partial contribution to its full size adds no complexity. Suppose the two subgraphs are 
cascaded in reverse order, the complexity becomes Q + QP. But Q + QP * P + PQ if P * Q. 
Thus, the order of the subgraphs is relevant to the complexity. If P<Q, P + PQ is always 
smaller than Q + QP. Therefore, the order selection should be such that the first one 
requires less complexity than the second one. The numbers of multiplications required for 
different combinations of u and v are shown in Table 2 and Table 3 for 4 x 4 and 8x8, 
respectively. Note that for u=0 or M/2 and v = 0 or N/2, no multiplication is required 
(except a left shift operation by one bit). 

Table 2. The Number of Multiplications Table 3. The Number of Multiplications 
Associated with Fuv (4 x 4) Associated with Fuv (8 x 8) 


X 
u 1 

0 

1 

2 

3 

0 

0 

3 

0 

3 

1 

3 

6 

3 

6 

2 

0 

3 

0 

3 

3 

3 

6 

3 

6 


\y o 

u \ 

1 

2 

3 

4 

5 

6 

7 

0 

0 

5 

3 

5 

0 

5 

3 

5 

i 

5 

20 

10 

20 

5 

20 

10 

20 

2 

3 

10 

6 

10 

3 

10 

6 

10 

3 

5 

20 

10 

20 

5 

20 

10 

20 

4 

0 

5 

3 

5 

0 

5 

3 

5 

5 

5 

20 

10 

20 

5 

20 

10 

20 

6 

3 

10 

6 

10 

3 

10 

6 

10 

7 

5 

20 

10 

20 

5 

20 

10 

20 
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According to Table 3, we can estimate the average number of multiplication requ red for 
a nonzero element in the 8 x 8 case. Assume that the chance of a nonzero element falling 
in any u-v pair is equally likely. Then the average will be the 1/64 of the sum of all the 
numbers shown in Table 3. The result is 9.5 multiplications per nonzero element. Th\s is 
about 7 to 20 times faster than the naive approach mentioned earlier. Similarly, from a 
2, we will get 3 multiplications per nonzero element for the 4 x 4 case. 

The update of a goal matrix is straightforward because it involves only additions of the 
corresponding elements in each partial contribution. The total number of additions is (n- 
1)MN for the update of the goal matrix, where n is the number of nonzero elements in an 

M x N target matrix. 

5. Performance Comparison of the Approaches 

The advantages of the proposed approach are as follows: 

(\) it has essentially no delay time and computational redundancy. 

(2) It allows any scanning or transmission patterns, including the zig-zag scanning pattern. 
Note that the zig-zag scanning pattern is generally good for many images. However, a 
better or optimal scanning pattern for a particular image may deviate from the zig-zag 
scanning pattern [14]. Furthermore, it can be different from one image to another. 
Therefore, it is critical to be adaptive to different scanning patterns. 

(3) In inverse PIT, it has lower computational complexity than traditional fast algorithms. 

Both advantages 1 and 2 are due to the separate processing of the element in the input 
matrix. The performance comparison of approach 1, approach 2, and the proposed approac 
are summarized in Table 4. 


i t • * tv* 



IHiBM 

Scanning Pattern 

| 

Fixed and Poor 

1 

Delay Time 


Low 

Lowest 

£ 

Computational Redundancy 

High 

Low 

Lowest 

XL. -* 

# of Multiplication* 

384 

114 

95 

# of Additions* 

1864 

740 

576 


'The fast algorithm in [8] is used for approaches 1 and 2. The number of nonzero elements 
in an 8 x 8 target matrix is assumed to be 10 on the average. In addition, 4 stages of image 
reconstruction are assumed. 


6. Conclusion 

This paper presents a new and promising solution to the problem of heavy computation of 
IDCTs in inverse PIT process. When approach 2 was shown to have poor and discontinuous 
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image build up in 1990 [10], the research in fast progressive reconstruction for transform 
domain PIT schemes seems to be hopeless. With the proposed approach, both the fast 
progressive reconstruction and the pleasant image build up can be achieved simultaneously. 
This is an encouraging result for the research of transform domain fast progressive 
reconstruction. 
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